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Abstract 

The connection between the neutrinoless double beta ((/3/3)oi'-) decay effective Majorana 
mass, |Mee|, and the branching ratios of the decays — > l^l'^, 1,1' = e,fi, of the dou- 
bly charged Higgs boson is analysed within the Higgs Triplet Model of neutrino mass 
generation. We work in the version of the model with explicit breaking of the total lepton 
charge conservation, in which l^l'^ , I J' = e,fj,,T, are the dominant decay modes 
of H^^. It is assumed also that are relatively light so that they can be produced at 
LHC and the branching ratios of interest measured. Taking into account the current and 
prospective uncertainties in the values of the neutrino mixing parameters most relevant for 
the problem studied - the atmospheric neutrino mixing angle 623 and the CHOOZ angle ^13, 
and allowing the lightest neutrino mass and the CP violating Dirac and Majorana phases to 
vary in the intervals [0,0.3 eV] and [0, 27r], respectively, we derive the regions of values of 
BR(^r±± ^ e±e±) and BR(i?±± e±A*^) for which |Mee| > 0.05 eV, or \Mee\ < 0.05 eV. 
This is done for neutrino mass spectrum with normal ordering, inverted ordering and in the 
case when the type of the spectrum is not known, and i) without using the possible additional 
data on BK{H^^ ii^fi^), ii) using prospective data on BK{H^^ n^fj,^). In the latter 
case results for several values of BR(_ff^^ — > are presented. 
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1 Introduction 



Determining the origin of neutrino masses and mixing is one of the major challenges of future 
research in neutrino physics. It is well known that the existence of nonzero neutrino masses can 
be related to the presence of more complicated Higgs sector in the Standard Theory, involving 
additional Higgs fields beyond the single doublet field. Actually, it was realised a long time 
ago [H El [3] that a Majorana mass term for the left-handed (LH) flavour neutrino fields can be 
generated by SU(2)l x U{1)y invariant Yukawa couplings of two lepton doublet fields to a Higgs 
triplet field, carrying two units of the weak hyperchange, |y| = 2. Such a Higgs field has an 
electrically neutral, singly charged and doubly charged components. The Majorana mass term for 
the active flavour neutrinos arises when the neutral component of the Higgs triplet field acquires 
a nonzero vacuum expectation value (vev), breaking the SU{2)l x U{1)y symmetry. There are 
several possible realisations of this scenario. The realisation in which the global U{1)l symmetry 
associated with the conservation of the total lepton charge L is broken only spontaneously by the 
Higgs triplet vev [4J, was ruled out by the LEP data on the invisible decay width of the Z'^-boson. 
If, however, the U{1)l symmetry is broken explicitly in a manner that leads to a nonzero vacuum 
expectation value of the neutral component of the Higgs triplet fleld (see, e.g. 13 [6]) one obtains 
a viable model of neutrino mass generation. This model has been investigated in detail recently 
[71 El [9] and was shown to have a rich and physically interesting phenomenology owing to the fact 
that i) the couplings of the doubly and singly charged Higgs fields to the flavour neutrinos and 
charged leptons are proportional to the elements of the Majorana mass matrix of the (flavour) 
neutrinos, Min, and can be relatively large, and that ii) the physical doubly charged and singly 
charged Higgs flelds, and H^, can have masses in the range from ~100 GeV to ~1 TeV 

and thus can, in principle, be produced and observed at LHC. Point i) implies that the indicated 
couplings are determined essentially by the elements of the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix [TU] and by the neutrino masses. In [3 El E] it was shown that by 
studying the decays — > I'^l'^ , I, V = e, /x, r, it might be possible to obtain information on the 
absolute neutrino mass scale, on the type of neutrino mass spectrum (which can be, e.g. normal 
hierarchical (NH), inverted hierarchical (IH) and quasi-degenerate (QD)), and on the Majorana 
CP violating phases [11] present in the neutrino mixing matrix. 

In the present article we investigate the possibility to use the information on the neutrino mass 
spectrum and the Majorana CP violating phases from the measurements of the — > l^l'^ 
decay branching ratios, BR(i^^^ — > l^l'^), 1,1' = e,/x, in order to obtain predictions for the 
effective Majorana mass in neutrinoless double beta ((/?/3)o!/-) decay, |Mee| (see, e.g. p2]). Our 
study is motivated by the fact that most probably the searches for the doubly charged scalars 
and the decays — > l^l'^ will be carried out at LHC before the next generation of (/?/?) oz/-decay 
experiments will be operative. Among the different decay channels — > l^l'^, 1,1' = e, /x,r, 
the easier to observe are those with two electrons (positrons), two muons (antimuons), or an 
electron (positron) and a muon (antimuon), e^e^, /^^/^^ and e^/i^, in the flnal state. If the mass 
of does not exceed approximately 400 GeV, the branching ratios of the decays into 
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e^e^, e^//^ and /^^//^ can be measured at LHC with a few percent error [T3]. We will show that if 
the doubly charged Higgs bosons will be discovered at LHC and at least the three branching 
ratios BR(i7±± e^e^), BR(ff±± e^^^^) and BR(i7=t± fJ-^f^^) will be measured with 
a sufficient accuracy, one can obtain unique information on the (/3/3)o!y-decay effective Majorana 
mass |Mee|- This information will be extremely important, in particular, for the upcoming next 
generation of (/?/9)oi/-decay experiments. 



2 The Higgs Triplet Model 

In the Higgs Triplet Model (HTM) a I = 1,Y = 2 complex SU{2)l triplet of Higgs scalar 

fields is added to the Standard Model (SM) Lagrangian. In the 2x2 representation, the Higgs 
triplet field has the form: 

^ A" -A+/V2 ) ^ ' 

where A*^, and are neutral, singly charged and doubly charged scalar fields. In the 

flavour basis in which the charged lepton mass matrix is diagonal we are going to use throughout 
this article, a Majorana mass term for the LH flavour neutrinos can be generated (without the 
introduction of a right-handed neutrino fields) by the SU{2)l x U{1)y gauge invariant Yukawa 
interaction: 

C = hn i,J,L C (irs) A i^iL + h.c. (2) 

Here hi'i = hw, I', I = e,fi,T, are complex Yukawa couplings forming a symmetric matrix h, C 
is the charge conjugation matrix, T2 is a Pauli matrix for SU{2)l indices, and tpfi^ = {uil II)'^, 
I = e, fi, T, is the LH lepton doublet field. A non-zero triplet vacuum expectation value, (A'^) = 
7^ 0, gives rise to a Majorana mass matrix M for the LH flavour neutrino fields i^il: 

Mi,i = 2hi,i{A^) =V2hinVA, l',l = e,fi,T. (3) 

The requisite f a 7^ arises from the minimisation of the most general SU{2) x C/(l)y invariant 
Higgs potential [11 [6]: 



V = ($t$) + Ai ($t$)2 ^ ^2 iv(AtA) + As [Tr(AtA)]2 + A3 Det(AtA) 
+A4 ($^$)Tr(AtA) + A5 ($t^.$)Tr(Atri A) + (^-^^ (^^iraA^^) + /i.cj 



(4) 



= ((/)"'" (f)^)'^ being the SM Higgs doublet field. In eq. dH), > is the common mass of the 
triplet scalars. The choice < ensures that (c/)'^) = v/^/2 0, which breaks spontaneously 
SU{2) X U{1)y to U{1)q. In the model of Gelmini-Roncadelh [i] the term /i($^ir2Al"<I>) is absent, 
which leads for < to a spontaneous breaking of the global U{1)l symmetry associated 
with the conservation of the total lepton number. The resulting Higgs spectrum contains a 
massless Goldstone boson - the triplet scalar Majoron, J, and another light scalar, H^. The 
decay — > H^J would give too large a contribution to the invisible decay width of the Z'^-boson 
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and this model was excluded by the LEP data. The inclusion of the term /i($-^ir2A^<I>) explicitly 
breaks the lepton number conservation when A is assigned two units of the total lepton charge L, 
and therefore avoids the presence of a Goldstone boson - the Majoron, in the model [21 E]- Thus, 
the scalar potential in eq. ^ together with the triplet Yukawa interaction of eq. ^ lead to a 
phenomenologically viable model of neutrino mass generation. 

The expression for v/\ resulting from the minimisation of the potential V, eq. (jlj), reads: 

''^A — TTTTn 7T , , o ) for IIA < f . (5) 

2Ml + (A4 + As) ' ^ ^ 

In the scenario of relatively light triplet scalars within the discovery reach of the LHC we will be 
interested in, one has Ma ~ v and eq. ([5]) leads to va ~ ^- In extensions of the HTM, the term 
/i($"^ir2A^$) can arise in various ways: i) through the vev of a Higgs singlet field [14^ I15j: ii) 
can be generated at higher orders in perturbation theory [6]; or iii) can appear in the context of 
theories with extra dimensions [5l [16| . 



An upper limit on v/\ can be obtained from considering its effect on the parameter p = 
M^/M^ cos^ Ow. In the SM, p = 1 at tree-level, while in the HTM one has 

p=l + 6p= , X = va/v. (6) 

The measurement p ^ 1 leads to the bound va/v < 0.03, or < 8 GeV. At the 1-loop level 
VA must be renormalised and explicit analyses lead to bounds on its magnitude similar to those 
derived from the tree- level analysis |17j . 

The HTM has seven physical Higgs scalar particles {H++,H — , H+ , , H° , A° ,h°). The 
doubly charged Higgs field H^^ coincides with the triplet scalar field A"*""*". The remaining Higss 
mass-eigenstate fields are in general mixtures of the doublet and triplet fields. The corresponding 
mixing parameter is proportional to the ratio of triplet and doublet vevs, va/v, and hence is 
small even if va assumes its largest value of a few GeV. Therefore ^ H^, are predominantly 
composed of the triplet fields, while is predominantly composed of the doublet field and plays 
the role of the SM Higgs boson. The masses of H^^, H^, H^, are of order Ma with splittings 
of order A5 v. For Ma < 1 TeV of interest for direct searches for the Higgs bosons at the LHC, the 
couplings hi'i are constrained to be 0(0. 1) or less by a variety of processes such as ^ — > eee, r — > III 
etc. These constraints are reviewed in \18\ I19j. 

In this article we will be interested, in particular, in the decays of into a pair of same-sign 
charged leptons, — > I'^l^, I', I = e, p, r, which give clear signals even in hadron colliders like 
the LHC. These decays are important not only for the searches of H^^, but also because of the 
very interesting possibility that their lepton flavour dependence can be directly related to the 
Majorana mass matrix of the LH flavour neutrinos. 

In our analysis we will assume that Mjj±± < Mjj± and f a ^ 1 MeV. Under these conditions 
the decay — > H^W^ is forbidden, while the decay — > W^W^ is sufficiently strongly 
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suppressed. In this case the branching ratios of the decays — > I'^l^ are given by the following 
simple expressions (see, e.g. OEllH]): 

BR,n^BR{H^^^l'^l^) = ^^^f^ (7) 

J- + Oi'i l^iq \rii'ir 

- (8) 



1 + ^I'l m: 



where 6i'i is the Kronecker delta. Note that the branching ratios depend only on the parameters of 
neutrino mass matrix. The measurement of BR;/; can give significant information on the elements 
of the neutrino mass matrix |M;/;|, and therefore, e.g. on the absolute neutrino mass scale (i.e. 
lightest neutrino mass), type of neutrino mass spectrum, Majorana CP violating phases in the 
neutrino mixing matrix, etc. 



3 The Neutrino Masses, Mixing and the (/3/5)oi/-Decay 

We work in the flavour basis in which the mass matrix of the charged leptons is diagonal. As we 
have shown, in the Higgs triplet model of interest, the LH flavour neutrino fields vil acquire a 
Majorana mass term. The corresponding Majorana mass matrix M is diagonalised with the help 
of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mixing matrix |1U| : 

Mw = [f7pMNsdiag(mi,m2,m3) C/pMNs]^', (9) 

where nij, j = 1,2,3, are the real positive eigenvalues of M//' - the masses of the Majorana 
neutrinos Xj with definite mass. In what follows we will use the standard parametrisation of the 
PMNS matrix (see, e.g. [201 [2l]): 



^ C12C13 S12C13 siae \ 

t^PMNS = 



-S12C23 - Cl2S23Sl3e*'' C12C23 - Sl2S23'Sl3e*'^ S23C13 



■ ^21 ■ 0^31 

diag(l, e*~, e* 2 ) 



\ S12S23 - Cl2C23Sl3e*'' -C12S23 - Sl2C23Sl3e*'' C23C13 J 



(10) 

where Cij = cosOij, Sij = smOij, the angles 9ij = [0,7r/2] (i < j = 1,2,3), 5 = [0, 27r] is the Dirac 
CP-violating phase, and a2i and asi are two Majorana CP-violation phases [111 I22j . The phases 
a2i and a^i can vary in the interval [0, 27r]. It proves useful for our further discussion to define 
also the difference of the two Majorana phases: 032 = 031 — 021- Let us add that at present we 
do not have experimental information on 5, Q21 and Q31. 

The existing neutrino oscillation data [23[|24[ 1^1261127) allow to determine with a rather good 
precision the mixing angles and neutrino mass squared differences which drive the solar neutrino 
and the dominant atmospheric neutrino oscillations, sin^20i2, Am2i and sin^2023) l^^-iiK— 
|Am32|), and to obtain a rather stringent limit on the CHOOZ angle ^13. In our analysis we will 
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use the following best fit values of sin 2(9i2, Am|i, sin^ 26^23 and lAm^^l [251 [2^150]: 

Amii = 7.6 X 10"^ eV^ sin^ 26^12 = 0.87, (11) 
|Am|i| = 2.4 X 10"^ eV^ sin^ 26^23 = 1 . (12) 

The upper limit on sin^ 20i3 obtained in CHOOZ reactor anti-neutrino experiment [26j reads: 

sin^ 26*13 < 0.14. (13) 

From the global analyses of the neutrino oscillation data one finds (see, e.g. [30]): 

sin^ %< 0.056 , 99.73% C.L. (14) 

The next generation of experiments with reactor I/g, which are under preparation, Dooble CHOOZ 
[31j . Daya Bay [32], RENO [33], can improve the currently reached sensitivity to the value of 
sin^ 2^13 by a factor of (5-10) (see, e.g. [M]); while future long baseline experiments aim at 
measuring values of sin 26*13 as smah as IQ-^-lO'^ (see, e.g. [35]1. 

Let us note that the uncertainty in the experimental determination of sin^ 2^23 corresponds to 
a rather large interval of allowed values of S03 [24J: 0.38 < s^g < 0.62. We will take into account 
this uncertainty in our numerical analysis □. It should be added that the accuracy on sin^ 2^23 
is planned to be improved considerably in future long baseline experiments. The uncertainty in 
sin^ 2^23 is foreseen to be reduced in the T2K experiment [36], for instance, to sin^ 2^23 > 0.99 
(0.45 < s^g < 0.55), if the true value of sin^ 2^23 = 1- As we will see, the correlations between 
the branching ratios of the decays — > l^l'^, BR^/^, 1,1' = e,/x, r, and the effective Majorana 
mass in neutrinoless double beta {{(3(3)ou-) decay, = M^e, which is the main subject of our 

study, depend not only on the elements of the neutrino mixing matrix, but also on the type of 
neutrino mass spectrum and on the absolute scale of neutrino masses. 

As is well known, owing to the fact that the sign of Arri^^, cannot be determined from the 
existing data, there are two possible types of neutrino mass spectrum compatible with the data - 
with normal ordering and with inverted ordering. In the standardly used convention we are also 
going to employ, the two spectra correspond to: 

- mi < 7712 < iTT'S, Arn-|]^ > 0, nor7nal orderi7ig (NO), 

— rns < mi < m2, ^rrisi < 0, inverted ordering (10). 

The v-mass spectrum can be: i) normal hierarchical (NH), mi <C m-2 < m^, with m2 — \J Amgi — 
8.8 X lO^'^ eV, m3 = y^^Ii — 4.9 x 10^^ eV; ii) inverted hierarchical (IH), m^ <^mi < m2, with 
7712 = a/AthI^ = 4.9 X 10~^ eV; mi = A/Am|^"^^"Am|^ = 4.8 x 10~^ eV, and iii) quasi- degenerate 
(QD), mi = m2 = m^, m\ 23^ |A?7i|]^|. In the latter case one has: mj ^ 0.10 eV. 

^Varying sin^ ^12 in the Scr interval of allowed values of sin^ ^12 [28l[30|, 0.25 ;$ sin^ S12 0.37, has essentially 
negligible effect on the results of our analysis. 
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The type of neutrino mass spectrum, i.e. sgn(Am|^), can be determined by studying oscil- 
lations of neutrinos and antineutrinos, say, <-> and i/^ <-> i>e, in which matter effects are 
sufficiently large. This can be done in long base-line i/-oscillation experiments (see, e.g. j35j). If 
sin^ 26*13 ~ 0.05 and sin^ 6*23 ^ 0.50, information on sgn(Am3]^) might be obtained in atmospheric 
neutrino experiments by investigating the effects of the subdominant transitions t'^(e) ^e(/i) 
and P^(e) ~^ ^e{n) of atmospheric neutrinos which traverse the Earth [3^. For v^[e) {or ^^i{e)) 
crossing the Earth core, new type of resonance-like enhancement of the indicated transitions 
takes place due to the (Earth) mantle-core constructive interference effect (neutrino oscillation 
length resonance (NOLR)) [38j [§. For Arrig^^ > 0, the neutrino transitions I'p.^e) ~^ ^e{fi) 
enhanced, while for Am|^ < the enhancement of antineutrino transitions ^'^(e) — > ^e(^) takes 
place, which might allow to determine sgn(Am|^). If sin^ ^13 is sufficiently large, the sign of Am|^ 
can also be determined by studying the oscillations of reactor z^g on distances of ~ (20 — 40) km 
[42j . An experiment with reactor v'g, which, in particular, might have the capabilities to measure 
sgn(Am|j^), was proposed recently in [33] (see also [Sj)- Information on the type of neutrino mass 
spectrum can also be obtained in /3-decay experiments having a sensitivity to neutrino masses 
~ ^J\Am^\ = 5 X 10~^ eV [35] (i.e. by a factor of ~ 4 better than the planned sensitivity of the 
KATRIN experiment [36], see below). 

Direct information on the absolute neutrino mass scale can be derived in '^H /3-decay exper- 
iments [47t I38| I46j . The most stringent upper bounds on the z^e mass were obtained in the 
Troitzk [48j and Mainz [36] experiments: 

mp^ < 2.3 eV, 95% C.L. (15) 

We have mp^ = mi^2,3 in the case of the QD z/-mass spectrum. The KATRIN experiment |46j . 
which is under preparation, is planned to reach a sensitivity of rrip^ ~ 0.20 eV, i.e. it will probe 
the region of the QD spectrum. 

The CMB data of the WMAP experiment [35] , combined with data from large scale structure 
surveys (2dFGRS, SDSS), lead to the following upper limit on the sum of neutrino masses (see, 
e.g. [50]): 

^ = S < (0.4-1.7) eV , 95% C.L. (16) 

j 

Data on weak lensing of galaxies, combined with data from the WMAP and PLANCK experi- 
ments, may allow S to be determined with an uncertainty of ~ 0.04 eV \50\ I51j. 

In our analysis we will consider both types of neutrino mass spectrum - with normal and with 
inverted ordering, as well as the specific cases of normal hierarchical (NH), inverted hierarchical 
(IH) and quasi-degenerate (QD) spectra. Correspondingly, the lightest neutrino mass min(mj) = 

^As a consequence of this effect, the corresponding i^^(e) {or J^^Ce)) transition probabilities can be maximal [31] 
(for the precise conditions of the mantle-core (NOLR) enhancement see [381 [39]). Let us note that the Earth 
mantle-core (NOLR) enhancement of neutrino transitions differs from the MSW one. The conditions of the 
Earth mantle-core enhancement [38 | l39 j also differ |40j from the conditions of the parametric resonance enhancement 
of the neutrino transitions discussed in the articles [41] . 
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Figure 1: Diagrams for the neutrinoless double beta decay. The diagram (a) is the standard and 
dominant one, and (b)-(f) are possible but negligible in the HTM. 



mo, which determines the absolute neutrino mass scale, will be varied in the interval: 

< mg < 0.3 eV , mo = min(mj) , j = 1,2,3. (17) 

As we will show, the results we obtain essentially do not depend on the maximal value mo as long 
as the latter is not smaller than 0.3 eV. The reason is that for mo ^ 0.3 eV (i.e. in the QD 
region), the branching ratios BR/// we are interested in practically do not depend on the neutrino 
masses: 

3(1 + Ovi) 

Neutrinoless Double Beta decay 

In the Higgs triplet model the massive neutrinos are Majorana particles. Determining the 
nature of massive neutrinos is of fundamental importance for understanding the origin of neutrino 
masses and, more generally, for understanding the symmetries governing the particle interactions. 
The existence of massive Majorana neutrinos is associated with non-conservation of the total 
lepton charge. In this case the neutrinoless double beta decay {A, Z) {A, Z + 2) + 2e~ is 
allowed (see, e.g. \12\ [52} I53j). Assuming that the dominant mechanism for the decay is the 
exchange of light Majorana neutrinos (Fig. HJa)), the half-life Ty^, for the decay is given by 

T^J^ = {G°''\M°''\^\Mee\^)~^ , (19) 

where G^"^ is a phase space factor and M^'^ is the nuclear matrix element of the process. All the 
dependence of T^'/2 on the neutrino masses and mixing parameters factorises into the (/3/3)o,y-decay 



Ui'jUij 



mo 



> 0.3 eV . 



(18) 
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effective Majorana mass M^e'- 



\Mee\ 



(20) 



The most stringent upper bound on |Mee|, \M(,e\ < (0.35- 1.05) eV was obtained by using the 
lower hmit T^/^ > 1.9 xlO^^ yr (90% C.L.) found@ in the Heidelberg- Moscow ™Ge experiment |55j . 



The IGEX collaboration has obtained the result T^Y^ > 1.57 x 10^^ yr (90% C.L.), from which 



-1/2 

the limit |Mee| < (0.33 — 1.35) eV was derived |5BJ. A positive (/3/3)oi/-decay signal at > 3a, 
corresponding to T^J^ = (0.69-4.18) x 10^^ yr (99.73% C.L.) and implying |Mee| = (0.1-0.9) eV, 
is claimed to be observed in [57], while a recent analysis reports evidence at Ga of (/3/3)o!/-decay 
with |Mee| = 0.32 ±0.03 eV |58j. Two experiments, NEM03 (with ^o^Mo, ^^Se, etc.) [59] and 
CUORICINO (with ^^Oxe) [60], designed to reach a sensitivity to \Mee\ ~ (0.2 - 0.3) eV, set 
the limits: \Mee\ < (0.61-1.26) eV p] and \Mee\ < (0.19 - 0.68) eV [60] (90% C.L.), where 
estimated uncertainties in the NME are accounted for. The two upper limits were derived from the 
experimental lower limits on the half-hfes of ^°°Mo and ^^^Te, T^Y^ > 5.8 x 10^^ yr (90%C.L.) [59] 



and T^Y2 ^ '^•^ ^ ^^'^'^ (90%C.L.) [60]. Most importantly, a large number of projects aim at a 



-1/2 

-1/2 

sensitivity to |Afee| ~ (0.01-0.05) eV [6l]: CUORE (^^Oxe), GERDA (^^Ge), SuperNEMO, EXO 
(i36xe), MAJORANA (^^Ge), MOON (^o^Mo), COBRA (^^^Cd), XMASS (^^^Xe), CANDLES 
(^^Ca), etc. These experiments, in particular, will test the positive result claimed in [57j . 



The predicted value of \Mee\ depends strongly on the type of i/— mass spectrum [S2l[2n]) more 
precisely, on the type of hierarchy neutrino masses obey. The existence of significant and robust 
lower bounds on \Mee\ in the cases of IH and QD spectra [62] (see also [63]), given respectively |l| by 
|Mee| ^ 0.01 eV and \Mee\ ^ 0.03 eV, which lie either partially (IH spectrum) or completely (QD 
spectrum) within the range of sensitivity of the next generation of (/3/3)ojy-decay experiments, 
is one of the most important features of the predictions of |Mee|. At the same time we have 
|Mee| ^ 5 X 10^^ eV in the case of NH spectrum |64j . The fact that max(|Mee|) in the case 
of NH spectrum is considerably smaller than min(|Mee|) for the IH and QD spectrum opens 
the possibility of obtaining information about the type of z^-mass spectrum from a measurement 
of |Mee| 7^ [Mj. More specifically, a positive result in the future generation of (/3/3)o^-decay 
experiments with \Mee\ > 0.01 eV would imply that the NH spectrum is strongly disfavoured (if 
not excluded). For Am|^ > 0, such a result would mean that the neutrino mass spectrum is with 
normal ordering, but is not hierarchical. If Am|^ < 0, the neutrino mass spectrum should be 
either IH or QD. 



^In the quoted upper bound for \Mee\ a factor of 3 uncertainty in the relevant NME (see, e.g. [54]) is taken into 
account. 

*Up to small corrections we have in the cases of two spectra \Mee\ ^ |Am|2 cos 20i2| (IH) and \Mee\ ^ 

mocos26i2 (QD). The possibility of cos20i2 = is ruled out at ~ 6a by the existing data which also imply that 
cos26li2 ^ 0.30 (0.26) at 2a {3a) We also have 2.07 x 10"^ eV^ ^ |Ami2| ^ 2.75 x 10"^ eV^ at 3a. 
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4 Prediction for |Mee| from Measurements of BR(iif 



±± 



In this Section we investigate within the HTM the predictions one can obtain for the (/3/3)o!y-decay 
effective Majorana mass |Mee| by using data on BR^/^. The dominant mechanism of (/3/3)oz/-decay 
- the Hght Majorana neutrino exchange, corresponds to the diagram in Fig. [TJ^a), the contributions 
from the diagrams in Figs. [T][^b)-(f) being negUgible [65^ 166]. 

We use three branching ratios, BRge, BRg^, and BR^^, in our analysis. The expressions for 
these branching ratios in terms of neutrino masses, neutrino mixing angles and CP violating 
phases read (see [3 El [9]): 



(E,m2)BRe 



(E.m2)BR 



■MM 



Cl2Cl3(-Sl2C23 - Ci2S23Sl3e jmi 
+ Sl2Ci3(ci2C23 - Si2S23Sl3e*'')m2e*°2l + S23Cl3Sl3m'3e*^"^' "''^ 
(-S12C23 - Cl2S23Sl3e*'^)^"T'l 



+ (C12C23 - Sl2S23Sl3e*'^)'m2e*"^l + 523^^3 "^36*°^^ 



2 „2 



(21) 



(22) 



(23) 



Given the solar and atmospheric neutrino oscillation parameters and the CHOOZ angle, 
sin^0i2, Am2i, sin^6'23, jAmg^^l and ^13, |Mee| depends on mo = min(mj), a2i5 ^31, b and on the 
type of neutrino mass spectrum (NO or 10). In the case of spectrum with 10 or of QD type, the 
dependence of |Mee| on 6*13 is relatively weak and can be neglected, as long as cos 2^12 ^ sin^ ^13. 
In this case |Mee| does not depend on the Majorana phase 031 and on the Dirac phase b. From 
the measurement of the three observables, BRge, BRe^, and BR^j^, three parameters, say, mo, 021 
and 031, can, in principle, be determined and information on the type of neutrino mass spectrum 
- with NH, IH or QD can be obtained [71E1I9]. This would allow to tightly constrain |Mee|. Let us 
review briefly the predictions for BRgg, BRg^, and BR^^ in the cases of NH, IH and QD spectrum 
(see also [HIHllS]). 



a) NH spectrum, mi ^ m2 < ms. 

Using the best fit values of the neutrino oscillation parameters one finds |§ that in this case BRee, 
BRe^i, and BR^^ can take values in the following intervals: ^ BR^^ ^ 10^^, ^ BRg^^ ^ 0.08, 
0.16 < m?^^ < 0.31. We get BR^^" = for (a32 - 25) = tt and sj^ = sj^V ^^21/ ^^'ii- 0-05, 
while BR^« = for (032 - 6) = ir and 5^3= sf2cf2 cot^ 6'23(Amii/Am|i)^ 6.9 x lO'^. The 

^The inequality cos 29i2 sin'^ ^13 is fulfilled for the 2a experimentally allowed ranges of values of cos 2^12 and 
sin^Sia, see. e.g. [281 130) . If one uses the 3a ranges, one obtains sin^ S13 / cos 2^12 ^ 0.22 . 

®The limiting values quoted in this paragraph are obtained for the best fit values of the neutrinos oscillation 
parameters and for sin^ 2^13 < 0.14. 
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minimal and maximal values of BRJJ^ depend weakly on sf^. Neglecting this dependence, one 
obtains a simple expression for the Majorana phase (difference) 032 in terms of BR 



1 

31 \ ^ -^^V/^ ~ ^23 



2 \Am21J C^2''23'^23 



cosa32 = ^ ^ -^P'r-^ (24) 



b) IH spectrum, 1x13 ^ mi < m2. 

We find very different results in this case: O.Scfg cos^ 26*12 ^ BR™ < O.Scfg or 0.06 < BR™ < 0.5, 
< BR™ < 0.48, < BR™ < 0.14. One has BR™ = for 5 = 0, 021 = tt and sfg ^ 0.036. Now 
BR™ exhibits a very weak dependence on s^g. For the Majorana phase Q21 we obtain in terms of 
BR^ 



ct - 2BR™ , , 

cosa2i = l- , (25) 

^''13''21*21 



c) QD spectrum, mi 2,3 ~ 0.1 eV. 

It is not difficult to convince oneself that the branching ratios of interest for the QD spectrum 
to a good approximation can take values in the following intervals: cos^ 2^12/8 ^ BRg° ^1/3 
or 0.03 < BR;^° < 0.33, < BR^° < 0.46, cos^ 26*23/3 < BR);^° < 0.33. Actually, we have up to 
small corrections BR|?g° = (2/3)BR™. For the Majorana phase Q21 in this case we get: 

cosa2i = l- 2 • (26) 

ZC]^3C2i621 

Given 021 and a sufficiently large S13, information about the Dirac phase 6 and the Majorana 
phase asi can be obtained from the knowledge of BRg° and BR^°. If, however, a stringent limit 
on S13 will be obtained, 031 can be determined using BR^° and the knowledge of 021- 

It is clear from the above simple analysis that the measurement of the three branching ratios 
BRge, BRg^, and BR^^ would provide information about the type of neutrino mass spectrum 
and the Majorana phases. If, for instance, it is experimentally established that BRee > 0.01, the 
neutrino mass spectrum of NH type will be excluded. The spectrum can either be of IH or QD 
type. If in addition BR^^ is determined to be BR^^ > 0.14, the IH spectrum will be ruled out. 
If, however, the neutrino mass spectrum will turn out to be QD, it will be very difficult (if not 
practically impossible) to distinguish between the spectrum with NO and that with 10, i.e. to 
get information about the sign of Am^i. 

Consider next the more general case of mg having an arbitrary value. First, let us set ^13 = 
for simplicity. We will consider the case of ^13 7^ later. For ^13 = 0, the main uncertainty in the 
prediction of | M^e \ comes from the lack of knowledge of mo and a2i . Note that in this case BRee 
and BRe^ are independent of 031, similarly to Mge Lj. Knowing these two branching ratios allows 

^BRer is independent of 031 as well, but this mode is more difficult to measure than the two modes we are 
discussing. 



10 



to determine niQ and cosa2i- In the case of spectrum with NO we have for the hghtest neutrino 
mass: 

mi Amai (2c?aCooBRee - s^oBReu) - m?BReu , , 

COSQ21 = , , \ '^P, \ ' . • (28) 



y^m^ + Arrigi 2c^2'5i2"ii\/"^i + Am2i(BRe^ + 2BReeC: 



Equation (I27p is vaHd also for the second to hghtest neutrino mass mi in the case of spectrum 
with 10. The expression for cos 021 obviously cannot be used to determine cos 021 for mi = 0: 
for ^13 = and mi = 0, neither BRee nor BRg^ depend on 021 • In the case of spectrum with 10 
(inverted ordering) we obtain: 

^2 ^ (2Ami3 - Am|i)(2ci3BRee + BR^^) - 2ci3(Am|3 - cf^Amj^) ^^^^ 



COSQ21 = 1 o-^r^ 9-T + O — T-] ■ (30) 

2cf2s22(BRe/, + 2BReec23) \l^mlj ^ ' 

As can be expected, for m\ ^ Am^i, the expression for cos 021 in the case of spectrum with NO 
coincides with that for spectrum with 10. Using eq. ()27p . we get a universal expression for |Mee| 
valid for both types of spectrum - with NO and 10 and any hierarchy between neutrino masses: 

,9 I sr^ 9 \ 2c9oAmoi — 2c9q(3s?9 — l)Am9i ^ , , 

\Mee\ =\ymj\ BRee = 9 9 ^BRee 31 

sgnfAmii ) X BR^p ■? 9 , ^ 

- 2.4 X 10"^eV^ (32) 

1 — 3BRee — 3BRe^ 

where in the last equation we have used the best fit value of ^23 and have neglected the term 
~ (35^2 — l)Am2i/Am|i. Note that the denominator in the expression for |Mee|, eq. (f3T|l . does 
not go through zero since we have: 

2ci3BRee + BRe„= 2ci3 ("^1 + ^12 Ami,) NO spectrum, (33) 
^■^ 3mf + Am|i + Amf, > > \ J 

24BRee + BRe^ = ^l^^ji^^^ , 10 spcctrum , (34) 
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where we have neglected terms ~ (Am2i/Am23) in the second equation. We see that in the 
QD region, where |Mee| has a relatively large value, one has (2ci3BRee + BRe^) ^ (2/3)4(1 + 
Ami3/(3m2)). 



It follows from eq. (|32|) that |Mee| > 0.05 eV ~ -y/jArnlJ can be predicted without the 
knowledge of sgn(Am3i), if the collider experiments show that the branching ratios BRee and 
BRe/i satisfy 

- ^BRee + ^ > BRe, > -^BRee + ^ (35) 
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If indeed sin^ ^13 is negligibly small and these conditions are satisfied by the measured BRee 
and BRe^, a positive result can be expected in the next generation of (/3/3)o;y-decay experiments 
having a sensitivity to |Mee| > 0.05 eV. Note that the magnitude of the left and right sides 
of the inequality is very sensitive to the value of C23. Note also that these conditions do not 
depend explicitly on BR^^. For this reason we will first obtain constraints on |Mee| using only 
the branching ratios BRgg and BRg^. 

Next, we analyse the case of 6*13 7^ numerically. We calculated BRee, BRg^, and Mee by using 
|Am3]^|, Am||, and sin^ 2^12 given in eq. (fTTI) and (fT2]) . We allow mo to vary in the interval in 
eq. (II 7|) . while the other parameters are varied in the following ranges reflecting the uncertainties 
in their knowledge or lack of any constraints: 

sin^ 26^23 > 0.94, sin^ 26*13 < 0.14, 5, 021, 031 = 0-27r. (36) 

Later we will present results for the prospective smaller uncertainties in sin 26*23 and sin^ 26*13, 
corresponding to sin^ 2^23 > 0.99 and sin^ 2^13 < 0.04. 

In Fig.[2]we show the regions in the BRge— BRe^ plane where we definitely have |Mee| > 0.05 eV 
or I Mee I < 0.05 eV. More specifically, the solid (red) line determines the complete allowed region 
in the HTM, corresponding to Arrigi, |Am|^| and sin^ 20i2 given in eqs. pi|) and p2|) . and values 
of mo, sin^2023) sin^20i3, 6, 021 and 031, which were allowed to vary in the ranges specified in 
eqs. (fT7|) and ([36]) . The dashed blue (dash-dotted green) lines determine the black (grey) regions 
where |Mee| is definitely larger (smaller) than 0.05 eV in the HTM when mo, sin^2023) sin^20i3, 
(5, a2i and a^i, are varied within the indicated intervals (i.e. eqs. (fT7|) and ([36|) ). For values of 
BRee and BRg^ from the region depicted in white (and located between those shown in black 
and in grey), the determination of |Mee| is not unambiguous: both values of |Mee| > 0.05 eV and 
I Mee I < 0.05 eV are possible. This degeneracy can be lifted to certain extent, but not completely, 
by using additional information on BR^^ (see further). The dotted black line in Fig.[2]corresponds 
to BRee + BRe^ = 1. We show it only to indicate the boundary of the region of possible values of 
BRee and BRe^: the area above this line is unphysical. The results in Fig. [2] are obtained without 
using the possible additional data on[§ BR^^. The left and middle panels correspond to NO and 
10 spectrum, respectively, while the results shown in the right panel were obtained assuming that 
the sgn(Am|;^) (i.e. the type of the neutrino mass spectrum) is unknown. The black area where 
I Mee I is, e.g. definitely larger than 0.05 eV in the right panel corresponds to the intersection of the 
black areas in left and middle panels. Note that we can have |Mee| ^ 0.05 eV also in the region 
shown in white and located between the grey areas in the right panel of Fig. [2j This cannot be 
unambiguously predicted, however, knowing only the values of BRee and BRe^^ which lie in the 
white area. 

Next we show how the results discussed above change when we add information on BR^^. 
The /i^^^ decay mode of is relatively easy to measure at LHC by virtue of the two same 
sign muons in the final state. We present results for BR^^ = 0; 0.1; 0.2; 0.3 in Figs. El HI [5] and 
m respectively, where BR^^ = in practice corresponds to BR^^ < 0.01. When we quote a 

*We recall that in this analysis we do not use possible data on BR^/j, BRt-t and BRer- 
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NO 10 unknown sign of 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 2: Values of BRge and BRg^, for which \Mee\ > 0.05 eV (black areas limited by the dashed 
blue lines) or |Mee| < 0.05 eV (grey areas limited by the dash-dotted green lines) in the HTM. 
The solid (red) line shows the entire region of allowed values of BRee and BRe^ in the HTM (black 
and grey areas and the white area between the coloured one). The results shown are obtained by 
varying ttt-q, sin^ 2^23) sin^ 2^13, 6, a2i and a^i, in the ranges given in eqs. (I17p and (I36p . The left 
and middle panels correspond to NO and 10 spectrum, respectively, while the right panel was 
obtained assuming that sgn(Am|^) is unknown. The dotted line corresponds to BRge + BRg^ = 1. 
The region above this line is unphysical. See text for further details. 



NO, BR^^=0 10, BR^^=0 unknown sign of Amg^, BR^^=0 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 3: The same as in Fig. [2l but assuming that the experimentally determined BR^^ = 0. 
The dotted line corresponds to BRee + BRe/^ + BR^^ = 1; the region above the line is unphysical. 
See text for further details. 
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Figure 5: The same as in Fig. [3l but for BR^^ = 0.2. 
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Figure 6: The same as in fig. [3l but for BR„„ = 0.3. 
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Figure 7: The same as in Fig. [21 but allowing sin^ 2^23 and sin^ 20i3 to vary in the following more 
narrow intervals: sin 2023 > 0.99 and sin^ 26*13 < 0.04. No information on BR^^ was used in 
deriving the results shown in the figure. See text for further details. 
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Figure 8: The same as in Fig. [71 but assuming that the experimentally determined BR^^ = 
The dotted line corresponds to BRee + BRg^ + BR„,j = 1; the region above the line is unphysical 
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Figure 9: The same as in Fig. [SI but for BR„„ = 0.2 
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Figure 10: The same as in fig. [H but for BRo„ = 0.3. 



specific value of BR^^ = x, we include an uncertainty of ±0.01 in x, i.e. we use BR^^ = j; it 0.01 in 
the numerical calculations. The dotted lines in Figs. [3E] correspond to BRge + BRg^ + BR^^ = 1. 
We do not use this constraint: the line represents the boundary of the physical region of values 
of BRee, BRe^ and BR^^. 

It is clear from Fig. [3] that the measurement of BR^^ can improve the predictability of |Mee|: 
the relative magnitude of the white "degeneracy" region, in general, is smaller than in the case 
when no information on BR^^ is available. This is not the case, however, for 10 spectrum and 
values of BR^^ = 0.1; 0.2 (Figs. |l]and[5l middle panels). 

Figures [2][6] show also the regions (dash-dotted line) where one definitely has |Mee| < 0.05 eV. 
If the measured values of BRee, BRe^ and BR^^ lie in one of these regions, the observation of 
(/3/3)o;/-decay in the next generation of experiments can be extremely challenging. Even in such 
a case, however, searches for the (/3/9)oy-decay are important and necessary also as a test of the 
HTM itself. If the (/3/9)o!/-decay is observed while the measured values of the leptonic decay 
branching ratios imply, e.g. a negative result of the searches for (/?/3)o,y-decay, we will be led to 
conclude that Mi'i and hin are not directly related: Mi'i ^ V^hi'iVA- Such a situation can arise, 
for instance, if va = 0, or in models with H~^^ which is not an SU{2)l triplet, but, e.g. is a y = 4 
singlet |67j with couplings to the charged leptons given by hin (I'j^)'^' Iff . 

We have performed the same analysis, but with reduced uncertainties in sin 26*23 and sin^ 26*13: 
sin^ 2^23 > 0.99 and sin^ 2^13 < 0.04. The indicated precisions (or better ones) in the determina- 
tion of sin^ 2023 and sin^ 2^13 are expected to be achieved in the upcoming T2K [36] and reactor 
antineutrino experiments Double CHOOZ [ST], Daya Bay [32] and RENO [33], respectively. The 
results of the analysis are shown graphically in Figs. ITlllOl The notations in Figs. [TlfTOl are the 
same as in Figs. [2M1 We see from Figs. [TlllOlthat improving the precision on sin^ 2023 and sin^ 2^13 
leads to a noticeable reduction of the regions of values of BRee and BRg^, for which it is impos- 
sible to determine whether |Mee| > 0.05 eV or |Mee| < 0.05 eV. As Fig. [H] demonstrates, the 
reduction will be particularly significant if the measur ed BR^^ < 10-2 (which we remind the 
reader corresponds to the case denoted by us as BR^^ = 0). 
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5 Conclusions 



We have investigated the connection between the (/3/3)oi/-decay effective Majorana mass |Mee|, 
and the branching ratios of the decays H^^ —>■ l^l'^, 1,1' = e, /x, of the doubly charged Higgs 
boson within the Higgs Triplet Model (HTM) of neutrino mass generation. Our analysis 

was performed within the version of the model with explicit breaking of the total lepton charge 
conservation, in which — > l^l'^, I, I' = e, /j, r, are the dominant decay modes of H^^. In this 
model the couplings of the doubly charged Higgs field H^~^ to the flavour neutrinos and charged 
leptons are proportional to the elements of the Majorana mass matrix of the (flavour) neutrinos, 
Min, and the branching ratios BR(i7^^ — > l^l'^) are entirely determined by the elements of 
the PMNS matrix and neutrino masses. The latter possibility is realised if the mass of the 
doubly charged Higgs scalar does not exceed the mass of the singly charged one, Mfj±± < M^± , 
and if the vacuum expectation value of the neutral component of the Higgs triplet field satisfies 
f^A ^ 1 MeV. The model under discussion was shown [71 El E] to have a rich and physically 
interesting phenomenology owing to the fact that the physical doubly charged and singly charged 
Higgs fields, and H^, can have masses in the range from ~100 GeV to ~1 TeV and thus 

can, in principle, be produced and observed at LHC. More importantly, by studying the decays 
e,^,T, it might be possible to obtain information on the absolute neutrino 
mass scale, on the type of neutrino mass spectrum and on the Majorana CP violating phases 
present in the neutrino mixing matrix [TJ El [9] . 

In the present article we have investigated the possibility to use the information on the neutrino 
mass spectrum and the Majorana CP violating phases from the measurements of the — > 
l^l'^ decay branching ratios, BR(//^^ — > l^l'^), 1,1' = e,id, in order to obtain predictions 
for the effective Majorana mass in neutrinoless double beta {{(3l3)ou-) decay, |Mee|. Among the 
different decay channels — > l^l'^, 1,1' = e,fi,T, the easier to observe and measure the 
corresponding branching ratios with high precision are those with two electrons (positrons), two 
muons (antimuons), or an electron (positron) and a muon (antimuon), e^e^, e^/i^ and /i^/i^, 
in the final state. If the mass of does not exceed approximately 400 GeV, the branching 

ratios of the decays into e^e^, e^/i^ and /^^/i^ can be measured at LHC with a few percent 
error |13j . 

Taking into account the current and prospective uncertainties in the values of the neutrino 
mixing parameters most relevant for the problem studied - the atmospheric neutrino mixing angle 
023 and the CHOOZ angle ^13, and allowing the lightest neutrino mass and the CP violating Dirac 
and Majorana phases to vary in the intervals [0,0.3 eV] and [0,27r], respectively, we have derived 
the regions of values of BR(ff^^ — > e^e^) and BR(i?^^ — > e^fj,^) for which we definitely have 
\Mee\ > 0.05 eV, or |Mee| < 0.05 eV. This is done for neutrino mass spectrum with normal 
ordering (NO), inverted ordering (10) and in the case when the type of the spectrum is not 
known. In what concerns the branching ratio BR(ff^^ — > fx^fj,^), we have considered two cases: 

i) the possible data on BR(i?^^ — > /i^//^) is not used as an additional constraint in the analysis, 

ii) the possible data on BK{H^^ fi^fj,^) is included in the analysis. In the latter case, results 
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for several values of BR(i?^^ /i^^u^) have been obtained. 

Our results are presented graphically in Figs. [2]|10[ They show that if the doubly charged Higgs 
bosons will be discovered at LHC and at least the two branching ratios BR(ff^^ — > e^e^) 
and BR(i7^^ — > e^fi^) will be measured with a sufficient accuracy, one can obtain important 
information on the (/5/3)oi/-decay effective Majorana mass \Mee\- In the various cases considered, 
we have identified the regions values of BR(i?=^± e^e^) and BR{H^^ e^^^), for which 
|Mee| is definitely bigger or smaller than 0.05 eV (Fig. [2]). We have shown also that due to i) the 
uncertainties in the determination of sin^ 2^23 and sin^ 2^13 , ii) the absence of data on the CP 
violating phases in the neutrino mixing matrix, and iii) the existing rather loose upper bound 
on the absolute neutrino mass scale, there exist also noticeable regions of values of BR(ff^^ — > 
e^e^) and BR(i?^^ — > e^/x^) for which it is impossible to determine unambiguously whether 
\Mee\ > 0.05 eV or \Mee\ < 0.05 eV (Fig. [5]). This "degeneracy" can be partially hfted by using 
the additional information from a measurement of BR^^ (Figs. [HE]). 

The same analysis was performed with reduced uncertainties in sin 26*23 and sin^ 26*13 cor- 
responding to sin^ 2^23 > 0.99 and sin^ 2^13 < 0.04. The results are presented graphically in 
Figs. [TlllOl They show that improving the precision on sin^ 2^23 and sin^ 2^13 leads to a notice- 
able reduction of the regions of values of BRee and BRg^ for which it is impossible to determine 
whether |Mee| > 0.05 eV or |Mee| < 0.05 eV. The reduction will be particularly significant if the 
measured BRuu < 10~^. 
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